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ABSTRACT 

The highly stable propagation characteristics of VLF 
transmissions make them an ideal source of world-wide 
navigation data. The arrival phase of a VLF signal is 
continuously compared with a local standard. The change in 
the phase difference between the received signal and the 
local standard is converted to relative motion with respect 
to the VLF transmitter. 

This investigation concerns itself with the development 
as a whole of the VLF navigation system. Refinements are 
made in the antenna system and the antenna signal processing. 
In addition, navigation procedures, tables, and forms are 
presented for system use. 
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I. 



INTRODUCTION 



Very low frequency transmissions form the basis of a 
very accurate long range navigation system. The character- 
istic of these transmissions that make them ideal for such 
a use is that their arrival phase at a point distant from 
the transmitter is both stable and predictable. 

The VLF wave propagates through a waveguide formed on 
one side by the ionosphere and on the other side by the 
surface of the earth. The phase of the received signal 
(the carrier frequency) at some point distant from the trans- 
mitter remains constant as long as the propagation path 
(waveguide) or the distance from the transmitter does not 
change. With the onset of darkness, the ionosphere height 
increases and modifies the waveguide, causing a change in the 
observed arrival phase. However, the observed arrival phase 
will return to its original value when the entire propagation 
path is again in daylight, and the waveguide assumes its 
initial dimensions. 

With the addition of two stable oscillators, one at the 
transmitter to serve as a reference for the master oscillator, 
and the second at some remote receiving site, motion of the 
receiver can be detected by observing the phase of the 
received signal as compared to the local oscillator. By 
commencing to track a transmitter signal in this manner from 
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a known geographic position, the navigator continuously has 
one line of position. Proper selection of a second trans- 
mitter will result in a second line of position crossing the 
first, and an electronic fix. 

The next section presents the background of this system, 
and the basic work accomplished in its development. Chapter 
III goes on to present an explanation of mode theory, which 
describes the propagation mechanism. This discussion leads 
to the basis of the system stability, and hence its accuracy 
Three basic problems are encountered in the use of this ' 
navigation system. These are examined in general in Chapter 
IV. The problems are long path signal contamination, pre- 
diction of the diurnal shift, and the conversion of phase 
information into geographic positions. Each is treated 
respectively in Chapters V, VI, and VII with their solutions 
Chapter VIII sets forth the system, and treats the units 
from an individual component view as well as a discussion of 
the system as a whole. Chapter IX presents the proposed 
procedures for use of the system. 
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II. BACKGROUND 



The development of the long range VLF navigation system 
starts with Blackband [1] and Stanbrough [2] , who both used 
this method experimentally with good results. Woods Hole 
Oceanographic Institute has developed a table of geographic 
positions [3] which contains bearings and distances to various 
VLF transmitters. While this table appears to be the initial 
effort in this line, the geographic positions are each one 
degree of latitude and longitude apart, a spacing too great 
for general navigation use. 

Lake [4] conducted a general investigation of this VLF 
navigation system, covering the basic propagation, phase- 
tracking receivers used, and general navigation procedures. 

McKay and Preston [5] covered in depth the factors which 
degrade the system accuracy, and developed the idea of using 
an antenna with a cardiod-shaped radiation pattern to reject 
the long path signal which could contaminate the desired 
signal in some cases. 

Roeder [6] developed an antenna system producing a 
cardiod-shaped pattern utilizing two loops and a whip with 
the associated equipment necessary to electronically rotate 
the resulting pattern. 

This paper will bring together the efforts of Stan- 
brough, Lake, McKay, Preston, and Roeder to result in a system 
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for long range VLF navigation. The specific results of this 
investigation are the development of a practical shipboard 
unit for use as a long range VLF navigation system, an 
improved method for predicting the time and amount of diurnal 
shift, the formulation of a means to generate navigational 
tables of positions, and the presentation of specific forms 
and procedures for system utilization. 
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III. MODE THEORY 



The propagation of electro-magnetic radiation in the 

VLF range is best described by waveguide mode theory. The 

guide, which channels or contains the wave, is formed by the 

earth on the lower side, and the ionosphere on the upper 

side, where the term "ionosphere" refers to the ionospheric 

_ 0 

height for which the conductivity is. equal to 10 mhos/m. 

As in conventional waveguides, many modes, other than 
the principal mode determined by the physical geometry, 
exist in the near region of sources and discontinuities. The 
resulting signal sensed at any point in the field is equal 
to the vector sum of the individual fields due to the multiple 
modes present. The higher order modes have greater attenu- 
ation factors than the principal mode, and therefore their 
contributions to the resultant signal decrease as the observer 
moves away from the transmitter. Table I gives various 
values of attenuation factors for modes 1, 2, and 3 of a 20 
kHz signal [7] . 



MODE ATTENUATION FACTOR 
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1.5 
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km 






night 
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km 
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night 
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km 
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night 


7.5 


dB/1000 


km 



TABLE I 

VLF ATTENUATION FACTORS 
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As an illustration, examine the 20 kHz signal at 1000 
km, 2000 km, and 4000 km from the source, during darkness. 

At 1000 km, the principal mode (n = 1) is down by 1.6 dB. 

The second mode is down by 2.7 dB, and the third mode is down 

by 7.5 dB. At 2000 km, the modes are down respectively 
3.2 dB, 5.4 dB, and 15 dB. At 4000 km the principal mode is 

down by 6.4 dB, the second mode is down 10.8 dB, and the 

third mode is down 30 dB. Further increases in distance 
reduce the signal contributions of the higher order modes 
and they may be neglected from consideration. In general, 
at distances greater than 4000 kilometers, the effects of 
higher order modes (n = 2 and greater) can be neglected. 

VLF propagation has also been described in terms of ray 
theory. Both mode theory and ray theory will provide the 
same results, however the mode theory is better suited in 
this case. 
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IV. BASIC PROBLEMS 



There are three major problems to be solved before the 
system can be used for navigation. They are the prediction 
of the time of occurrence and amount of diurnal shift, the 
means of nulling or canceling the long path signal when the 
receiver is in the multipath zone, and the problem of conver- 
ting phase difference information into geographic motion 
and position. Each of these areas is discussed briefly below, 
and examined in detail in following sections. 

A. DIURNAL SHIFT PREDICTION 

Consider an example where the receiving site is located 
to the west of the transmitter and both places are in day- 
light. The recorded phase difference between the received 
signal and a local oscillator should have a constant value 
since there is no motion toward or away from the transmitter, 
and the propagation path is not changing. With the occurrence 
of sunset at the transmitter, the recorded phase difference 
at the receiver starts to increase as the height dimension 
of the earth-ionosphere waveguide increases. The increasing 
height causes a decrease in the phase velocity. The decreased 
velocity leads directly to the change or increase in the 
observed phase difference. The phase difference continues 
to increase until sunset occurs at the receiving site. This 
increase in the phase difference is the diurnal shift. With 
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the propagation path completely in darkness, the phase 
difference is constant until sunrise occurs at the trans- 
mitter. This time the phase difference decreases until it 
reaches its daylight value where it again tracks constant. 
Since long range navigation at sea is a 24 hour requirement, 
it is highly desireable that the amount and time of the 
diurnal shift be predictable. This in fact can be achieved 
and will be discussed in a later section. 

B. LONG PATH SIGNAL 

Utilizing the virtual ionosphere-earth waveguide and the 
low attenuation factor for the principal mode, the signal 
radiated from the antenna has three prominent characteristic: 
regions, each being associated with certain sections of the 
"guide." In an area (circle of radius 4000 kilometers) 
around the transmitter, the multi-mode zone exists. This 
region is not suitable for phase tracking due to the many 
modes present. Therefore, this area is not useable for 
navigating. From the multi-mode zone outwards from the 
transmitter is the single mode zone which is useful for phase 
tracking, and likewise for navigation. The single mode zone 
extends to a region which surrounds the antipode of the 
transmitter. Herein is the multi-path zone. This region 
exists since the transmitter radiates its energy omni- 
directionally. Therefore, in the area of the antipode, two 
signals are possible, the direct or short path signal, and 
the long path signal, that is one that has travelled more 
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than halfway around the world. The presence of these two 
signals can contribute major navigational errors if they are 
not processed with care. The situation can go from a mild 
contamination of the desired signal to complete dominance 
of the undesired signal. Of particular concern is the case 
where minor errors could be introduced into the receivers 
by the presence of this long path signal, and these would 
not be immediately detected by the navigator since the system 
would continue to track in the proper direction. 

The degree of contamination will be determined by the 
lengths of the two signal paths, and the signal attenuation 
factors for each path. In general, seawater paths have 
smaller attenuation factors than land paths, and daylight 
paths have greater attenuation than darkness paths. There- 
fore, a combination of darkness and seawater for the long 
path can provide a signal of comparable strength to that of 
the short path where it may have a significant portion over 
land and in daylight. 

It is highly desirable to use the multi-path zone in that 
it is also a single mode zone, and has very stable signals 
for phase tracking. If this region was not used, the system 
would be significantly restricted. The use of an antenna 
system with a directive pattern which discriminates against 
the undesired signal would solve this problem. 

C. TABLES 

In order to convert the phase change (receiver information) 
into geographic motion and position, suitable tables, similar 
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to those existing for LORAN A, must be generated to allow 
plotting the ship's position on charts. These tables should 
include known geographic positions, and bearings and 
distances (expressed both in nautical miles and microseconds) 
to transmitting stations. The geographic positions should 
be spaced no greater than thirty miles apart. 

D. MINOR PROBLEMS 

Four minor problems associated with maintaining an 
accurate plot of a ship's position are: 

a. Cycle slipping 

b. Transmitter outage 

c. Frequency offset of the transmitter 
master oscillator 

d. Sudden ionospheric disturbance (SID) . 

Each of the above should pose no problem for a navigator 
with a minimum familiarity with this system. Cycle slipping 
refers to the situation when the receiver is tracking a 
changing arrival phase condition, and drops behind in 
indicated phase by 360 degrees or integer multiples thereof. 
This would generally occur only where the phase difference 
is changing rapidly. This could only happen at the onset 
or offset of the diurnal shift. The type of receivers used 
in this investigation allows the operator to "drive" the 
receiver back onto the proper cycle if slipping does occur. 
The operator also has an option to choose a shorter time 
constant for the receiver to insure that cycle slipping 
does not occur. 
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Scheduled transmitter outages are no problem. Normally, 
these outages are scheduled, and the navigator should have 
shifted the receiver to an alternate transmitter. Tracking 
of a transmitter can commence from any known position. By 
selecting an optimum time, the shift of stations can be 
made without difficulty, and without a loss in position 
accuracy. Since a prudent navigator values a three lines of 
position fix, risk from this and other anomolous outages 
can be minimized by tracking three transmitters at all times. 

Frequency offsets of transmitter master oscillators ■ are 
scheduled sufficiently in advance that the amount of the 
offset and date and time of its application will be promul- 
gated to all concerned. The navigator at sea need only 
make an appropriate correction for his calculations. 

Sudden ionospheric disturbances are caused by solar flares, 
magnetic storms, and high altitude nuclear blasts. The 
duration and magnitude of the distrubance will depend on the 
duration and magnitude of the disturbing phenomena. In 
general the disturbances are not intense, and of short length. 
In any case, their effect is to lower the ionospheric height, 
causing an increase in the phase velocity of propagation. 

This increase would make the transmitters being tracked 
appear closer to the ship. Continuous chart recorders would 
help the navigator to recognize the onset of a distrubance . 
There would be a rapid decrease in the displayed phase 
difference information. Under this circumstance, a dead 
reckoning track should be plotted from the last accurate 
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position. Subsequent positions may be determined by 
celestial navigation or other means . These positions 
compared with the system performance will provide an indica- 
tion of when the ionospheric behavior has returned to normal. 

E. PROPAGATION CONSIDERATIONS 

The navigator, in using this system, must understand 
certain VLF propagation characteristics to intelligently 
assess the validity of signals under various circumstances. 
For example, it must be realized that the velocity of signal 
propagation is dependent on the height of the ionosphere 
and also on the signal path. The signal path is classified 
as to being over seawater, land, or ice. If the path is over 
ice, the ground conductivity is such that the wave behaves 
as if it were propagating over land. Hereafter, reference 
will be made only to the land path due to the similar 
behavior of propagation over ice. Since this is a relative 
navigation system, the percent of signal path over land or 
seawater has no effect on the system ability to fix a 
position during daylight hours. However, as it will be seen 
in a following section, the observed diurnal shift for the 
land portion of the path is not proportional to that for the 
seawater path. Therefore, the navigator must be cognizant 
of the amount of propagation path over seawater, and over 
land to accurately predict the diurnal shift corrections. 

If a polar path is being utilized, he must also be knowledge- 
able of the extent of the ice pack for that particular time 
of the year. 



18 



As an example of the differences in velocities of signal 
propagation, TABLE II demonstrates their magnitudes for a 20 
kHz signal under various conditions. The four starred 
conditions denote the standard ionospheric conditions. Using 
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PHASE VELOCITIES UNDER VARIOUS CONDITIONS 

these velocities one arrives at the standard conversion value 
of 6.18 microseconds per nautical mile for daytime propagation 
over seawater. Night propagation over seawater has a conver- 
sion of 6.22 microseconds per nautical mile. 

Navigation tables are calculated later in this investiga- 
tion. The entries are on the basis of 6.18 microseconds per 
nautical mile. They may be adjusted for the second conversion 
(6.22) by multiplying the distances expressed in microseconds 
by a factor of 1.006. 
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V. ANTENNA SYSTEM 



The ability to discriminate against the long path signal 
when the ship is located in the multi-path zone is the 
function of the directive antenna system. An ideal system 
should have small, lightweight, and stationary elements. 
However, normal multi-element arrays utilized as antennas 
in the VLF range would be of such size that only a land-based 
installation could support such configurations. 

The use of a loop antenna, whose vertical plane is pointed 
at the transmitter, provides a means of attaining the above 
favorable antenna properties with bi-directional sensitivity. 
When this pattern is coupled with a properly phased whip 
signal, a cardiod-shaped pattern results with a null opposite 
to the bearing of maximum signal gain. 

The loop antenna output is the difference between the 
signals detected by opposite legs of the loop. A signal or 
voltage difference exists due to the spatial separation 
of the two antenna legs. If the antenna's orientation is 
changed, the magnitude of the detected difference signal 
decreases as the spatial separation (and hence the phase 
difference) with respect to the wave front decreases. 
Therefore, a null occurs when the plane of the antenna is 
perpendicular to the wave front. 

To attain the cardiod shaped antenna radiation (field 
intensity) pattern, the whip antenna signal is matched in 
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magnitude to that of the loop antenna, but with a 90 degree 
phase lag. The result of adding the two antenna outputs is 
a signal of twice the strength in the preferred direction, 
and a null 180 degrees from the maximum. Figure 1 is a 
composite of the radiation patterns of the single loop antenna, 
the whip antenna, and the resulting cardiod shaped pattern. 

Since the cardiod-shaped pattern is the result of a phase 
addition, and the navigation system utilizes changes in 
phase difference to detect motion, the change in phase across 
the axis of the pattern is of major concern. McKay and 
Preston [5] demonstrated both theoretically and experimentally 
that a "zero phase change window" does exist on the pattern 
axis. This investigation also confirms the presence of this 
window. Since this window is of finite width, errors will 
not be introduced in the system if the antenna pattern is 
not pointed exactly at the transmitter. If the resolver 
rotor indicates a direction within 5 degrees of the transmitter 
bearing, no appreciable error exists. McKay and Preston also 
have shown that the pattern null is down 30 dB for a sector 
of 40 degrees width. This investigation again confirms this 
performance. Therefore, any observed phase changes are 
indicative of the motion of the receiver toward or away from 
the transmitter and not due to a phase change introduced by 
the antenna system. 

The use of a single loop and single whip for phase track- 
ing is impractical in that only one station may be tracked 
for each antenna installation, and also a sizeable and 
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Figure 1. Composite Antenna Patterns. 
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massive antenna drive system would be required to continuously 
point the loop antenna. Roeder [6] developed a means by which 
2 stationary loops duplicate the performance of a rotatable 
single loop. As seen in Figure 2, the 2 loops (mounted 
concentrically and perpendicular to each other) drive field 
coils mounted in quadrature. The field generated in each coil 
is proportional to the signal sensed by the corresponding 
antenna. The rotor will sense the maximum signal when 
positioned at an angle that corresponds to the direction of 
signal arrival. The antenna radiation pattern illustrated 
in Figure 2.c is obtained by turning the rotor through 360 
degrees. This antenna system duplicates the behavior of a 
single loop pointed directly at the transmitter. 

As a demonstration of the feasibility of this approach, 
Roeder used goniometers in place of the field and sense coils 
to test a system capable of tracking two transmitters. A 
means of positioning the goniometer rotors was required to 
keep the antenna patterns pointed at the transmitters while 
the ship maneuvered. This was accomplished by use of a gyro 
compass signal input to a power amplifier which fed a servo 
motor attached to the rotor shafts. With the addition of a 
whip antenna and cardiod units, the proper radiation patterns 
were formed. 

In part, this investigation was directed toward the 
realization of a more practical shipboard installation. To 
reach this aim, the six feet diameter loop antennas used by 
Roeder were replaced by a submarine VLF crossed loop antenna. 
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Figure 2. a. Basic Schematic. 




Figure 2.b. Geographic Figure 2.c. Antenna 

Picture Pattern. 

Figure 2. Effective Single Loop Antenna System. 
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This greatly reduces the size of the topside or mast-mounted 
unit. An order of magnitude decrease is possible since the 
submarine antenna loops are wound on ferrite cores. 

In place of goniometers, precision resolvers were used 
with a reverse signal flow. Again, as in the case of the 
antenna, a much smaller physical size is realizable for an 
operational unit, and as a result allows the direct use of a 
gyro repeater signal for driving or positioning the rotors 
without the need of a power amplifier. 

The gyro compass repeater input can be either 120 VAC, 

60Hz or 120 VAC, 400 Hz. The synchro receiver drives a local 
compass card. Two other compass cards are geared to the 
driven unit. A means of securing the rotors of the resolvers 
to the compass cards is provided through a bearing pointer. 

The bearing of the station to be tracked is indicated by the 
pointer when the unit is in operation. 

Audio amplifiers are used to boost the VLF signal level 
for the two loops' inputs into the resolvers. The two 
amplifiers feed four windings. Each amplifier feed corre- 
sponding windings, one in each resolver. Separate gain 
controls are present for each winding in each resolver. More 
resolvers may be fed to allow a capability to simultaneously 
track more transmitters. Caution should be exercised to 
insure that sufficient signal energy is available at the 
output of the amplifiers for this purpose. 

In order to combine the loop and whip signals properly, 
commercial cardiod units are used in this system. These 
units are manufactured by TRACOR, Inc. of Austin, Texas. 
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They provide a means to vary independently the magnitude and 
phase of the whip signal. The loop signal may be reversed 
in phase. The output of the two cardiod units furnish the 
input or antenna signals for the VLF tracking receivers. 

Figures 3, 4, and 5 represent the results of tests 
performed on the antenna system using VLF transmitters NAA, 

NBA, and NPM as radiating sources. Each of the figures 
contain two plots. A circular plot of the field intensity 
demonstrates in each case that the axis of the pattern points 
in the direction of the transmitter. Also, the cardiod-shaped 
patterns are well defined, and clearly exhibit the null and 
severe attenuation in the reverse direction. The second plot 
is the measure of the phase distortion error. These plots 
are centered on the bearing of the transmitter and span a 
sector of 20 degrees. These measurements confirm the existance 
of the "zero phase change window" on the pattern axis, and 
the proper formation of the antenna pattern. 

The tests were conducted in the following manner. Two 
receivers were used to track a transmitter. One receiver used 
the antenna system as a signal source, and the second receiver 
used a whip antenna only. This second receiver served as a 
reference for comparison with the performance of the antenna 
system channel. The field intensity was measured every 5 
degrees of rotation of the resolver rotor, and plotted as 
shown in Figures 3. a, 4. a, and 5. a. For the phase distortion 
error measurements, the phase difference data of both receivers 
was compared as the resolver rotor position was changed. This 
information is presented for a 20 degree window in Figures 3.b, 
4.b, and 5.b. 
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Figure 3. a. Intensity Pattern. 
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Figure 3.b. Phase Distortion Error in 
Microseconds . 



Figure 3. Antenna Field Intensity Pattern and 
Phase Distortion Error Source - NBA, 
Bearing 115. 
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Figure 4. a. 
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Figure 4.b. Phase Distortion Error 
in Microseconds. 



Figure 4. Antenna Field Intensity Pattern 

and Phase Distortion Error Source - 
NAA, Bearing 060. 
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Figure 5.b. Phase Distortion Error 
in Microseconds. 



Figure 5. Antenna Field Intensity Pattern and 
Phase Distortion Error Source - NPM, 
Bearing 255. 
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VI. DIURNAL SHIFT 



The observed relative phase difference between the stable 
local oscillator and the received VLF signal (hereafter 
referred to as the "phase difference") at a remote receiver 
increases as the ionospheric height increases with the 
occurrence of sunset. The phase difference settles at a new 
constant value when the complete path from the transmitter 
to the receiver is in darkness. The reverse process of 
sunset is observed at sunrise. This change between the two 
constant values of the phase difference is called the diurnal 
shift. In order that the navigator may determine his position 
from sunset to sunrise, the diurnal shift must be predictable 
in the time of occurrence, and the amount of additional phase 
difference . 

McKay and Preston [5] approached this problem by calcu- 
lating various relative phase velocities for different 
conditions. A more direct solution to this problem is a 
combination of the work of Wait [8] and Brady and Crombie [9]. 
Wait's work, which is described in the next paragraph, gives 
us the amount of diurnal shift, and Brady and Crombie provide 
the means to calculate the time of sunset and sunrise at the 
height of the ionosphere. The work of Brady and Crombie will 
be treated more fully in a following section on navigation 
procedures . 

This solution assumes that the daily phase difference 
variation for a stationary receiving site can be accurately 
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represented by a trapezoid as shown in Figure 6 . For this 




Figure 6. Diurnal Phase Shift. 



example, consider the case where the receiving site is to the 
west of the transmitter. During the period represented by 
the interval A to B, there is a constant phase difference 
observed. At the time B, sunset occurs at the transmitter, 
and the height of the ionosphere starts to increase causing 
an increase in the phase difference due to the decreasing 
phase velocity. The phase difference continues to increase 
until sunset occurs at the receiver. At this time, the 
waveguide path has again stabilized, and a constant phase 
difference is observed until sunrise occurs at the transmitter 
at time D. Then the phase difference decreases until time E 
when sunrise occurs at the receiver. The phase difference 
is again at the same value as for the preceding day. If the 
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receiver was to the east of the transmitter, point B would 
represent the time of sunset at the receiver and C, the time 
of sunset at the transmitter. Likewise, D would be the time 
of sunrise at the receiver, and E would be the time of sunrise 
at the transmitter. However, the measured value of At would 
be the same in both cases. 

Wait developed the following equation 

at = 2_ <d/v n - c/v a ) 

where At is the total diurnal phase shift in microseconds, 

D is the distance in kilometers between the transmitter and 
receiver, c is the velocity of light, v^ is the phase velocity 
of the signal at night, and v^ is the phase velocity during 
the day, the last three quantities having units of meters per 
second. 

The foregoing presupposed a stationary receiving site. 

Now, with the receiver located on a ship that is moving away 
from the transmitter, Figure 6 must be modified by adding 
point by point a second plot which is a line with a slope 
equal to N x 6.18 microseconds per hour where N is equal to 
the speed of the ship in knots in the direction of the 
transmitter. (Note, the slope is positive for motion away 
from the transmitter, and negative for motion toward the 
transmitter.) At the same time, the night or darkness portion 
of the trapezoid pattern will be modified since At varies 
linearly with D. Therefore, for a more complete picture of 
the diurnal shift as it applies basically to the VLF navigation 
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system. Figure 7 is presented with the various portions of 
the composite process. 




Figure 7. Diurnal Shift with Ship Motion. 



Examining Wait's equation, the distance to the transmitter 
must be known, and ratios of the velocity of light to the 
phase velocity for both day and night must be determined. 

The great circle distance can be calculated using methods 
such as Bowditch [10]. In order to find the velocity ratios, 
an equation from the Omega System handbook [11] is used: 
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v h . 

-E = 1 - 0.36 — 
c a 



(27m -V - V ) £ 

^ 4tt/ 2 f h. 



Where 

v 

P 

c 

0.36 



h. 

l 



signal phase velocity 

free space velocity of light 

constant chosen to give a good fit 

ionospheric height for v/hich the conductivity 
,-6 



n 

f 



is equal to 10 mhos/m 

radius of the earth = 6378.155 kilometers 
phase shift on reflection at the ground 
phase shift on reflection at the ionosphere 
mode number (n = 1 for the principal mode) 
frequency in Hz 
Initially, an attempt was made to compare the results 
of the observed diurnal shift with the calculated value of 
At using the above equations. Since the distance was fixed, 
it was possible to arrive at accurate values for (c/v - 
c/v d ) • The measured values did not agree with those 
calculated . 

A second attempt to predict the diurnal shift began with 
measuring At for stations where the propagation path was 
entirely over land (stations NAA, NSS , and NBA). This resulted 

The heights used 
for the ionosphere were 70 kilometers for day, and 90 kilo- 
meters for night. The constant, 0.36, was replaced by "C" , 
and the equation solved for "C". This resulted in three 
values for "C" as shown in Figure 8. As an initial interpre- 
tation, "C" was represented by a linear function of frequency. 



in three accurate values for (c/v “ c/v^) 
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Therefore, the constant for a good fit has been replaced by 
a linear function of frequency. Using the data from Figure 
8, a curve showing the value of versus frequency is 
plotted as Figure 9 for the land propagation route, where 
is defined as (0 . 53961/0 . 3) (c/v - c/v^) . This substitution 
allows one to use nautical miles directly for the computation 
of At. 

The same procedure was followed for sea propagation 
routes. These results are shown in Figures 10 and 11. A 
similar curve is shown, Figure 12, which was derived from 
data presented by Spies and Wait [7] . Figure 13 is an overlay 
of Figures 9, 11, and 12 to compare the general shapes and 
slopes of the three curves. 

Figures 9 and 11 along with a modified version of Wait's 
equation, would be used by the navigator in his prediction 
of 

At = D C 

v 

the amount of the diurnal shift. In the case where the entire 
propagation path is over seawater, C v is read from Figure 11 
for the proper frequency. Then, At would be equal to the 
product of C v and D. If the propagation path for the signal 
is over both land and seawater. At is then the sum of the 
diurnal shift for the land portion of the path, and the 
seawater portion of the path. For example, examine the 
situation where a ship is 6000 miles from the transmitter 
(f = 22.0 kHz) . The seawater portion of the path is 4000 
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Figure 12. Reference Curve 
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Figure 13. Comparison Curves 



nautical miles. The diurnal shift due to the land portion 
is equal to 2000 x 0.0133 or 26.6 microseconds. The diurnal 
shift due to the seawater portion is equal to 4000 x 0.0196 
or 78.4 microseconds. The total diurnal shift for this path 
is 105.0 microseconds. 
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VII. TABLES OF POSITIONS 



In order to use this system of navigation, tables of 
positions must be generated that will allow the navigator to 
convert the measured changes in signal phase to geographic 
coordinates . 

Since this system is based on relative changes in signal 
phase, some standard procedure must be adopted for the 
navigational tables. This investigation used the great 
circle distance in nautical miles from a known position to 
the transmitter. (This distance is also listed in micro- 
seconds to facilitate the plotting of new positions. The 
generated tables also list the bearing of the transmitter 
and a second bearing, 90 degrees from the first to aid the 
navigator in laying down lines of position on charts and 
plotting sheets.) 

The known points are spaced every fifteen minutes of 
latitude and longitude to reduce the need for interpolation 
over large distances. 

Contained in the Computer Output section is a sample page 
from the tables. Across the top of the page is the name of 
the transmitting station, its frequency, and geographic 
location. Latitude is given in the left most and right most 
columns, and longitude is denoted in the first row. For any 
given position, there are four entires. The first line is 
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the great circle distance in nautical miles to the transmitter. 
The second line contains the same distance expressed in micro- 
seconds. The third line contains two entries, the first is 
the bearing to the station along the great circle path, and 
the second is a bearing to be used to lay down a very good 
approximation to the circular line of position. The navigator 
will be concerned with the bearing of the transmitter along 
the great circle path, in order to have the "zero phase change 
window" of the cardiod shaped radiation pattern pointed in 
the proper direction. Also, when the ship is located in the 
multi-path zone, knowledge of this bearing will allow the 
radiation pattern null to be pointed at the long path or 
undesired path signal. 

The method of calculation of the tables is a FORTRAN IV 
program, contained in the Computer Program section which was 
run on an XDS 9300 computer located at the Naval Postgraduate 
School. The mathematical model was derived from Ref. [10]. 

The program computes two quantities, the distance and true 
bearing from a known geographic position to a transmitter. 

The calculations are based on the following two equations. 

hav D = hav DLo cos cos l^ + hav 1 

hav C = sec esc D (hav col^ - hav(D~coL^)) 

The function hav (haversine) of an angle is equal to the 

2 1 

sine squared of half the angle (hav A = sin 2 A) . The quan- 
tities in the equations are as indicated below. Figure 14 
shows the basic navigational triangle as it appears in this 
case . 
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Figure 14. Solid Geometry 

Navigational Triangle. 
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D 



DLo - 



L 

L 

1 



1 

2 



C 



CoL 2" 



D~coL^- 



distance expressed in degrees between the known 
position and the transmitter 

difference in longitude measured from the known 

position to the transmitter, and identified as 

either east or west (This quantity will be less 

than 180 degrees.). 

latitude of the known position 

latitude of the transmitter 

difference in latitude (1 is measured from the 
known position latitude to the transmitter 
latitude and identified as either north or south.), 
an angle indicating the direction of the trans- 
mitter from the known position labeled north or 
south to agree with and east or west to agree 
with DLo 

colatitude of the transmitter (C 0 L 2 = 90 - are 
of the same name, and C 0 L 2 = 90 + L 2 if L^ and 
L 2 are of contrary names) 

the numerical difference between D and coL^ 



The computer program processes the angular position data 
in two formats, degrees and radians. The radian format is 
necessary for use in the computer angular subroutines. The 
distance is first calculated and then used as an input for 
arriving at the bearing to the transmitter. The program 
includes the necessary logic for determining the correct 
procedures and labels in the computation of DLo, 1, C 0 L 2 , plus 
the conversion of C into a true bearing. 
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For the purposes of this investigation, tables were - 
generated for an area of the Pacific Ocean adjacent to the 
coast of California, centered around Monterey. (The coordinates 
of this area are from 34N to 40N, and from the coast out to 
122W.) These tables are in the same format as that in the 
Computer Output section, and are in the custody of the Dean 
of Research Administration at the Naval Postgraduate School, 
Monterey, California. Similar tables may be produced for 
any other area of interest by using the proper geographic 
locations as input data for the computer program. 
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VIII. VLF NAVIGATION SYSTEM 



Herein is a listing of the navigation system components 
and a discussion of their properties and characteristics 
as they apply to the system. The units are referred to by 
their manufacturers and model numbers. These particular 
components were used in this investigation since they were 
available at the Naval Postgraduate School. Similar units 
from other manufacturers meet the requirements for this 
system and will provide the needed performance. 

A. COMPONENTS 

The basic shipboard system consists of the following items 

One stable oscillator 

Three VLF phase tracking receivers 

One VLF whip (sense) antenna 

One VLF crossed loop antenna 

One antenna control unit 

1. Stable Oscillator 

The very heart of the shipboard system is the stable 
oscillator. The stability of this oscillator will provide 
the mark and measure of the system accuracy. The transmitters 
listed in Table III all have a frequency stability of 1 part 
in 10^ (The meaning of "1 part in 10^" is that for a given 
oscillator or standard, the output frequency will deviate at 
the most 1 hertz in 10 11 hertz from the nominal value) [12] . 
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As an illustration of system accuracy, let us assume the 
shipboard oscillator has a stability also of 1 part in 10^. 
Now given the worst instance where these errors are additive, 

2 parts in 10^\ examine the system accuracy at this point. 
Since there are 0.864 x 10^ microseconds in one day, the 
above error of 2 parts corresponds to 2.32 microseconds. 

Using the conversion of 6.18 microseconds per nautical mile, 
and 2000 yards per nautical mile, the maximum error of 2 parts 
corresponds to 750 yards. If the shipboard oscillator has a 
stability of 2 parts in lO^, the corresponding maximum' 
error increases to 1125 yards. 

Present day frequency standards fall into two general 
categories, either a high precision temperature controlled 
crystal, or a device utilizing a natural atomic resonance. 
There are two stable oscillators located at the Naval Post- 
graduate School, one each of the types mentioned above, both 
of which are suitable for system use. 

The high quality quartz crystal oscillator is a Model 
2.5, manufactured by the Sulzer Laboratories of Rockville, 
Maryland. This unit requires a 22-32 VDC, 0.33 ampere input, 
with outputs (1 volt to 50 ohms) at frequencies of 2.5 mHz, 

1 mHz, and 100 kHz. The crystal unit operates at 2.5 mHz and 
is kept at a constant temperature by means of a two stage 
oven [13]. The crystal will drift with age, and monitoring 
of its output is necessary to determine the amount of drift. 

Prior to use of this type of crystal oscillator, it is 
necessary that the unit be in operation for periods from 3 
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to 6 months for settling. Loss of power to the unit will 
require the long stabilization procedure. Therefore, this 
type of oscillator is generally powered from wet cell batter- 
ies, which are constantly on charge from a normal power supply. 
On loss of power, the oscillator is then completely carried 
by the batteries . 

Atomic resonance type standards utilize the emitted 
energy of a given natural decay transition as a basis of 
frequency definition. The Model V-4700 Rubidium vapor fre- 
quency standard manufactured by Varian Associates of Palo Alto, 
California, uses a hyperfine transition of Rubidium 87 to 
stabilize a 5 mHz crystal oscillator [14] . This unit has 
outputs of 5 mHz, 1 mHz, and 100 kHz. It has a long term 
stability of 5 parts in lO"*"^, and a short term stability of 
1 part in 10'*''*'. The Varian model uses a similar power supply 
as the Sulzer crystal unit. The atomic resonance type 
standards require only a few days for the stabilization 
process as compared to the months required for crystal units. 

The Sulzer crystal has a short term stability of 8 
parts in 10'*''*', which would correspond to a maximum error of 
3375 yards. This amount of error is greater than an amount 
to be tolerated for the system. On this basis, a stability 
of 4 parts in 10'*''*' would provide accuracy of less than one 
nautical mile. 

In the selection of a frequency standard, one has the 
option of either type of the two discussed above or perhaps 
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a Cesium standard. Two performance characteristics must be 
evaluated, stability (or inversely dispersion) and drift. 

The Cesium standard has no drift but a moderate dispersion. 

The dispersion problem could be remedied by suitable averaging 
(integration) techniques, but the prime drawback for use in 
this navigation system is its high cost. Going to a Rubidium 
standard, the dispersion factor is improved (reduced), but 
drift now enters the picture. This adds the requirement for 
continual correction to maintain an accurate standard. Use 
of a highly precise crystal has the same problems as the 
Rubidium standard, but with greater drift. 

Therefore, with both of the practical standards, drift 
must be determined by comparison with a refined standard. 
Depending on the amount of drift and the smallest correction 
that may be made to the particular frequency standard, a 
regular correction procedure must be conducted to maintain 
an accurate navigation system. 

2 . VLF Tracking Receivers 

TRACOR Series 599 VLF Tracking Receivers, manufactured 
by TRACOR, Incorporated of Austin, Texas, were used in this 
investigation. These receivers were designed for time and 
frequency measurements with the VLF transmitting stations 
listed in Table III. 

Figure 15 is a simplified block diagram that depicts 
the operation of phase tracking VLF receivers. As illustrated, 
the output of the phase comparator is the input of the 
electronic phase shifter. The circuit is such that the phase 
error is driven to zero, and the system tracks an internal null. 
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Figure 15. Simplified Block Diagram of 
VLF Phase Tracking Receiver 



This receiver is powered either from 120 VAC (48-400 
Hz) or from a suitable DC source. Useful outputs of this 
receiver are as follows : 

a) time difference register - a digital counter that 
indicates the relative phase difference in microseconds 
between the received VLF carrier and the local frequency 
standard. 

b) audio - 1 kHz intermediate frequency. 

c) phase difference (100 microsecond scale) for display 
on an external chart recorder. 

d) relative carrier level from a front panel meter and 
also from a recorder output. 

The receiver requires a 0.01 microvolt signal at its terminals 
to enable normal phase tracking [15]. 

3 . Antennas 

The sensing apparatus for this navigation system 
consists of three separate antennas . Two loop antennas are 
concentrically mounted at right angles in a submarine VLF 
antenna. This antenna is designated as an AT-317C/BRR (fre- 
quency range of 14.6 - 38.0 kHz). For 'a shipboard installa- 
tion, one loop is aligned with the ship's centerline, and the 
other loop is at a right angle or running port and starboard. 
Each of these loops has a separate but identical cable run 
to the Antenna Control Unit. A VLF whip antenna (TRACOR Model 
No. 599800) is also used to provide a sense input to obtain 
the desired antenna radiation pattern. The only requirement 
placed on the location of the whip antenna relative to the 
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loop antennas, is that the distance should be small compared 
to a wavelength. This requirement provides no restriction 
when operabing at VLF frequencies. 

4 . Antenna Control Unit 

The Antenna Control Unit consists of the Antenna 
Control Console, the two cardiod units, plus the associated 
auxiliary components. Its purpose is to combine the signals 
from the three antennas into an effective single antenna 
signal. The resulting antenna pattern is to be rotatable to 
allow the ship freedom of motion to maneuver, and also to 
allow the navigator the freedom of choice to change the 
station being tracked. Therefore, a gyro compass repeater 
signal is utilized to provide the direction sense, and an 
antenna pattern control capability. An earlier section 
discussed the actual operation and performance of the 
components and the unit as whole. 

B. SYSTEM 

Figure 16 depicts the components of the navigation system 
and their interconnection. Two rotatable antenna radiation 
patterns (formed by means of the three antennas, the resolvers, 
and the cardiod units), sense VLF signals. These detected 
signals are compared in the phase-traking receivers with the 
local standard through the use of an internal frequency 
synthesizer and a phase-lock loop technique. The phase 
difference (a relative quantity) is read out on a digital 
indicator on the face of the receiver, and/or displayed on a 
continuous graph recorder as a function of time . The output 
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Antenna 
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Figure 16 . VLF Long Range Navigation System 

Components and Their Interconnections. 



format of the phase difference is in microseconds of change, 
which may be converted to distance traveled by the conversion 
unit of 6.18 microseconds per nautical mile, or through the 
use of tables. 

While Figure 16 shows RCVR #3 being fed from the whip 
antenna, a third resolver channel can be used to increase the 
system flexability. The whip signal is useful only in the 
single mode for a transmitter. 

The antenna radiation pattern must be capable of tracking 
the transmitter to insure that the received VLF signal enters 
through the zero phase change window on its axis. To achieve 
this capability, it is best to consider the components of the 
antenna patterns. The cardiod-shaped pattern is composed of 
the "figure eight" pattern from the "single effective loop" 
and the whip pattern. The cardiod orientation is the same 
as that of the loop antenna pattern, and this is controlled 
by the position of the resolver rotor. Therefore, the means 
of rotating the resultant cardiod-shaped pattern is to rotate 
the resolver rotor. 

As course changes are made, the ship "rotates" under the 
rotor of the gyro compass, which maintains its orientation. 
Likewise, the true bearing to a transmitter will not change 
with the ship's course change (the relative bearing measured 
from the ship's head or bow will change an amount equal to the 
course change). Therefore, by controlling the position of the 
resolver rotor with a gyro repeater input, the antenna pattern 
will remain pointed at the transmitter regardless of the 
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ship's maneuvers. Small corrections will be made as the ship 
traverses large distances, and the bearings to each trans- 
metter changes. Tables of known positions, discussed in the 
next section, will contain these bearings. The navigator need 
only make these changes to the resolver indicators. 
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IX. NAVIGATIONAL FORMS AND PROCEDURES 



A, FORMS 

Four new navigational forms are proposed for use with the 
VLF navigation system. Their purposes and use are explained 
in the following sections . 

Reference Position Form 

This form, Figure 17, records a position that is fixed 
by outside means. The form has a space to record the oscilla- 
tor drift in microseconds per day. The entry would be that 
determined by the last calibration check. Space is provided 
for a description of the method of fixing the ship's position, 
also. Typically, at the start of a voyage, an appropriate 
location can be used as a starting reference. At the time 
of the navigational fix, the receiver microsecond output is 
recorded. Space is provided for the recording of the geo- 
graphic line of position in microseconds from the table of 
positions. A is the numerical difference in the receiver 
output and the table of position value, such that A may be 
added to the receiver reading (also plus corrections) to 
obtain the proper chart value for plotting purposes. 

Hourly Navigation Record 

This form. Figure 18, provides the main log for hourly 
recording of the ship's location. The receiver output is 
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DATE 

TIME 

POSITION 



RCVR #1 RCVR #2 RCVR #3 

STATION STATION STATION 

RCVR RDG RCVR RDG RCVR RDG 

TABLE TABLE TABLE_ 

A A A 

OFFSET OFFSET OFFSET 

OSCILLATOR DRIFT 



DESCRIPTION OF FIX 



FIGURE 17. Reference Position Form. 
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entered in the column labelled "RCVR". The total correction 
to the receiver output is entered in the next column labelled 
"CORRT" . This should include A, oscillator drift, and At 
if appropriate for that particular time of the day. The hourly 
position may be entered in the two extreme right hand columns. 
Sufficient space is provided at the bottom of the form to 
allow the navigator to record the predicted values for At 
for night navigation. 

Sunset/Sunrise Calculation Forms 

These two forms. Figures 19 and 20, are used to calculate 
the time of sunset/sunrise at ionospheric heights . Both forms 
are self explanatory, and are based on the work of Brady and 
Crombie [9]. Appendix A is a collection of useful data for 
use with these two forms . 

B. PROCEDURES 

Navigation procedures can best be described in a narrative 
type discussion of a navigator's use of this sytem. 

Prior to the actual use of the navigation system, the 
navigator selects two or more transmitting stations for 
tracking. Among the things to be considered in his selection 
are bearings to the stations in the area where the ship will 
be steaming, scheduled transmitter maintenance periods, trans- 
mitter power, and distance to the station. Bearings and 
distance to the transmitting stations will be listed in the 
tables of position. Stations should be selected so that their, 
bearings will cross at least 30 degree angles (preferrable as 
close to 90 degrees as possible) to insure maximal accuracy. 
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E = -GHA 


E = 360 -GHA 




E = -GHA 


E 


= 360-GHA 




15 


15 




15 




15 




0 D = — 
U 15 






e D = A_ 

U 15 








+ 


— 




+ 




- 




12.0000 


xxxxxxx 




12.0000 




XXXXXXX 




E 






E 








a w 

6 d 


E 




s D " 




E 




LHA 


XXXXXXX 




LHA 




XXXXXXX 





TS TS 

(L and 6 are positive if North, and negative if South) 
FIGURE 19 . Sunset Calculation Form. 
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SUNRISE 



DATE 



L 






L 








A 






A 








6 






6 








GHA (1200) 






GHA (1200) 








SIN L 




A 


SIN L 






A 


COS L 




c 


COS L 






c 


SIN B 




B 


SIN 6 






B 


COS 6 




D 


COS 6 






D 


COS (LHA) =-.1276 


- A x B 




COS (LHA) =-.1276 




■ A X B 






C x D 






c 


X D 




LHA = (LHA) 
15 






LHA = (LHA) 
15 








GHA < 180 or 


GHA > 180 




GHA < 180 or 


GHA > 180 




E = -GHA 
15 


E = 360-GHA 
15 




E = -GHA 
15 


E 


= 360-GHA 
15 




0 D = - — 
U 15 






0_ = A 
D T5 








+ 


— 




+ 




— 




12.0000 


XXXXXXX 




12.Q000 




XXXXXXX 




E 






E 








w 

9 D 


E 




w 

9 D 




E 




LHA XXXXXXX 






LHA XXXXXXX 









TR TR 



(L and 6 are positive if North, and negative if South) 
FIGURE 20. Sunrise Calculation Form. 
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The distance to a transmitting station should be greater than 
2400 nautical miles to insure that the receiver is tracking 
the principal mode of the transmitting station. The Time 
Service Announcement of the U. S. Naval Observatory, a sample 
of which is reproduced as Table IV, is a necessary document 
for the navigator's use. 

It is also necessary that the navigator know the amount 
of frequency offset between his local standard and reference 
oscillators for the stations to be tracked. This offset 
should be expressed in microseconds and must be applied daily 
to the output of the receivers. As an exampled of this 
correction, suppose that the transmitter reference oscillator 
is ahead of the local standard by 1 part in 10^. This 
corresponds to a daily correction of 1.16 microseconds to be 
added to the output data (phase difference) of the particular 
receiver tracking that transmitter. Similar appropriate 
corrections will be applied to the other receiver outputs. 

After the stations have been chosen, the navigator will 
check the. antenna field intensity radiation patterns, and 
commence receiver tracking. As the ship gets underway, the 
necessary information is entered in the Reference Position 
Form . As time progresses, appropriate entires are made in the 
Hourly Navigation Record . As stated earlier the value of 
the oscillator drift was recorded on the Reference Position 
Form . Care must be taken to insure that this correction and 
that for frequency offset are properly applied. 



66 







/■—N. 














4-> 














CD 










U 




P 


£ 












VO 








CO 


CO 0 




X • 


U 


>i 












CD 






p 


u m 




N 




o u 


+> 


• 






iH 


i— 1 






O 


P 


0 


o pa 


P 


-p 0 


Cli Eh 














< — 1 *H 


o & 


co 


CD ^ 


0 


> 


0 


D 






U 


• 






rd co 


2 u 


rd 


• 




LD 0 


a 








<D 


o 






•H CO 




£1 r- o 


co 


LO 


X 


o 






o 


2 






U -pH 


O >1 ft H CO 


i — i 


0 


0 


o 






g 








CD g 


£ Xl 




• 


0 


g 4-> 










(D 








Qa CO 


4-> 


• 


P r- 


P 


0 P 


u 


CN 






4J 








CO P 


Td 


u 


O rH 


tr> 


U P 


p 








O. 








(d 


p 0 


p 




■H 


4H -H 


o 


O 




CP> 


(D 








u 


o £ 


o 


0 P 


co 


g JP 


-P 




VO 


CO 








Eh 


44 O p i — 1 O 




0 




















• — 1 




X 


0 


co .p 


p 


LO 




t — 1 


r- 










2 H 1 


0 


rd 4 -> 


g 


U 4-) 


0 


LO 






(N 










CO 0 


p 


-P O 


•pH 


O o 


> 


CO 




u 












Pm 44 


o 


co P 


Eh 


2 vo 


0 


CN 




(D 


44 


























X* 


0 


























0 






CO 














>1 








4-> 


• 




2 




(D 


rd 








rd 








U 


P 




O 




U 


o n3 






o 


TJ 








o 


fi 




M CO 




P 


o -H 






o 


P 










< 




co 2; 




(d 


cd Pi 






CD 


O 














CO o 




P 


rH Pm 






*H 


2 








H 1 


• 




H H 




(D 




















CO 




2 Eh 




4-> 


O 4P 






O Jp 














CO < 




P 


4-> O 






•P 


u 










• 




2 eh 




•H 


rd 








rd 






>4 O 




Eh 




< CO 




(d 


o 0 






o 


0 






Ph C* 








p< 




s 


o 






o 








O co 




CO 




Eh O 






^ £h 






CN 


Eh 






Eh O 




CD 




H 






«H D 






«H D 






<! CN 




r 0 




>h Q 




















> 




(D 




U < 




















• 




CO 




2 K 




^ £ 
















W O 




U 




W 




«H 0 4^ 


i — 1 






CN 








CO 




CD 




D P] 




(d-P^ 


V > 














CQ • 




a 




a 3 




P (d 
















O Q 




p 




w > 




•H -H P 


o 






O 












CO 




3 < 




g T3 0 


o 














^ - 




v — ’ 




Pm 2 




0 ® £ 


o 






H 








< 2 












2 O 
















> O 








Q • 




C4 


1 — 1 














C Eh 








2 co 


• • 


















2 O 








<3 


4-> 


















2 








• 


a 


















• H 








W D 


(D 


0 * 
















co pq 








2 


44 


p — 


o 






o 








co 






ro 


H 2 


44 


(D N 


CD 






o 








• < 








Eh O 


(D 


p pa 


• 






• 








D £ 






CO 


« 






r- 




















w 


Pm Pm 


P 


(D ^ 


i — i 






CN 














M 


O 


*H 


U 






















P$ 


Pm 




Pm 






















w 


W PJ 


CO 
























CO 


PI > 


*H 


























D 




























Q 2 


(D 




<r> 




















Eh 


W O 


i-4 




~ • 








• 












2 


w 


P 




VO 








CD 














u 






i — i 








ro 














CO 


CD 




0 o 








o 












W 




Jp 




P r- 








o 












o 




a 




*H VO 








r- 












2 




0) 




rd 








0 












P3 








2 - 








P - 












o 




Cn 


P 


2 








O 2 












2 




P 


O 


* cr> 






V. 


c<j ro 
















•H 




U -• 






rd 


^ • 












3J 




£ 


■P 


0 co 






O 


rH CO 
















0 


(d 


rH CO 






X 


rd o 












w 




1 — 1 


u 


4-> o 






rH 


P o 












o 




1 — 1 


o 


P ^ 






rd 


rd o 












H 




0 


PI 


U ^ 






PQ 


U o 












> 




44 
























p: 




























w 




(D 


p 






















CO 




Jp 


o 


























Eh 


•pH 






















pD 






-P 


< 






< 




















(d 


< 






PQ 














H 




• 


4-> 


2 






2 














Eh 




l — 1 


CO 

















-p 

0 P * 
4-> CO 43 • 

rd rd N 
2 o PP 

p 

0 • ^ ld 

P 0) CN O 
P 0 

•pH g P ^ 
-P *H 0 CN 
P 4-> 

O (DP 
U P H O 
(D 43 

2 ,P fd -P 
CO -P -P O 
h O U) P 



67 



Frequency Radiated Special 

Station Location (kHz)* Power (kw) Maintenance Transmissi 

NLK Jim Creek, Wash. 18.60 250 1000 to 1500 UT FSK continuou 
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TABLE IV - Time Service Announcement. 



The navigator will calculate the time of sunset for each 
of the transmitting stations using a Sunset Calculation Form . 
Using the ship's speed and course, he will estimate his 
position at ionospheric sunset, and calculate this time of 
sunset. These calculation will require the use of the Nautical 
Almanac to determine the declination and the 1200 Greenwich 
Hour Angle for the sun. Now using a modified version of 
Wait's equation for the amount of diurnal shift (the equation 
now allows the use of distance in nautical miles) we arrive 
at the value of At. 

At = D C 

v 

The value of D is readily accessible from the tables of posi- 
tion. Again using predicted course and speed, the various 
times of sunrise may be computed. C is obtained from Figures 
4 and/or 6, depending on the propagation path. Sufficient 
information is now at hand for the diurnal shift prediction 
trapezoid. The trapezoid must be modified to indicate a slope 
due to the motion of the ship with respect to the transmitting 
station. The amount of correction for each hour is N times 
6.18 microseconds, where N is the ship's speed component in 
knots toward (minus correction) or away from (plus correction) 
the transmitter. With the construction of the modified 
trapezoid completed, the hourly corrections are entered on 
the bottom of the Hourly Navigation Record for use at the 
required time. 

As an illustration of the accuracy of this method for the 
prediction of At, Table V lists the following results for two 
stations, NPM and NAA. 
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TABLE V 

TRAPEZOID PREDICTION RESULTS 



NPM 

frequency - 23.4 kHz 

propagation path - seawater 

predicted At - 41.2 microseconds 

mean error - 2.5 microseconds (810 yards) 

NAA 

frequency - 17.8 kHz 

propagation path - land 

predicted At - 30.3 microseconds 

mean error - 2.1 microseconds (680 yards) 

In the above two examples, mean error is equal to the average 
of the absolute values of the differences between the predicted 
trapezoid and the recorded diurnal shift. Samples of the 
differences were taken every fifteen minutes during the period 
that diurnal shift was present. 

Occasionally, it is necessary for transmitting stations 
to make offsets to their master oscillators to stay in tolerance 
with the national standard. These offsets are scheduled well 
in advance of their accomplishment. The navigator must make 
a similar correction to his readings at the same time as the 
actual offset is applied to maintain a tracking navigation 
system. This offset is applied only if navigation is in 
process. This correction can be ignored when starting from 
a new reference position, since in this case, the offset is 
then included in the A value. 
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APPENDIX A 



NPM 

L = 21-25. 5N 

X =158-09. 3W 
158.1550W 

A = .3652 
C = .9308 
0 D = 10.5437W 



NBA 

L = 09-03. 3N 

X = 79-38. 9W 
79.6483W 

A = .1574 
C = .9876 
© D = 5.3099W 



GBR 

L = 52-22N 

X = 1-11W 

1.1833W 

A = .7919 

C = .6106 

0 D = . 0789W 



STATIONS 



L 

X 

A 

C 



NSS 

L = 38-59 . IN 

X = 76-27. 2W 
76 .4533W 

A = .6291 
C = .7773 
0 D = 5.0969W 



NLK 

L = 48-12N 

X =121-55W 
121.9167W 

A = .7455 

C = .6665 

© D = 8.1278W 



NAA 

= 44-38. 9N 

= 67-16. 9W 
67.2817W 

= .7028 
= .7114 
= 4.4854W 



NWC 

L = 21-49 .OS 

X =114-09. 8E 
114.1633E 

A =-.3716 
C = .9284 
0 D = 7.6109E 



WWVL 

L = 40-40N 

X =105-02W 
105 . 0333W 

A = .6517 

C = .7585 

© -= 7.0022W 
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Time or X 
in Decimals 



Time or X 
in Minutes 



0000 


00 






0083 








0167 


01 






0250 








0333 


02 


Time or X 


Time or X 


0416 




in Decimals 


in Minutes 


0500 


03 


.0017 


.1 


0583 




.0033 


.2 


0667 


04 


.0050 


.3 


0750 




.0067 


.4 


0833 


05 


.0083 


.5 


0916 




.0100 


.6 


1000 


06 


.0117 


.7 


1083 




.0133 


. 8 


1167 


07 


.0150 


.9 


1250 ' 








1333 


08 






1416 








1500 


09 






1583 








1667 


10 






1750 








1833 


11 






1916 








2000 


12 . 






2083 








2167 


13 






2500 








2333 


14 






2416 








2500 


15 






2583 








2667 


16 






2750 
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Time or X 
in Decimals 

.2833 
.2916 
.3000 
.3083 
.3167 
.3250 
.3333 
. 3416 
.3500 
.3583 
.3667 
.3750 
.3833 
.3916 
.4000 
.4083 
.4167 
.4250 
.4333 
.4416 
.4500 
.4583 
.4667 
.4750 
.4833 
.4916 
.5000 
.5083 
.5167 
.5250 
.5333 
.5416 
.5500 
.5583 



Time or X 
in Minutes 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 
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Time or A 
in Decimals 

.5667 
.5750 
.5833 
.5916 
.6000 
.6083 
.6167 
.6250 
.6333 
.6416 
.6500 
.6583 
.6667 
.6750 
.6833 
.6916 
.7000 
. 7083 
.7167 
.7250 
.7333 
.7416 
.7500 
.7583 
.7667 
.7750 
.7833 
.7916 
.8000 
.8083 
. 8167 
. 8250 
.8333 
.8416 



Time or A 
in Minutes 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 
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Time or X 
in Decimals 

. 8500 
.8583 
.8667 
.8750 
.8833 
.8916 
.9000 
.9083 
.9167 
.9250 
.9333 
.9416 
.9500 
.9583 
.9667 
. 9750 ' 
.9833 
.9916 
1.0000 



Time or X 
in Minutes 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 



75 



£ 



NPM 23.4 KHZ 

LATITUDE 



NORTH 


114-00 


00-00 


2885.37 

17831.57 

299.4/029.4 


OC-15 


2878.03 

17786.22 

299.2/029.2 


00-30 


2870.74 

17741.15 

298.0/029.0 


00-4 5 


2863.49 

17696.36 

298.8/028.8 


01-00 


2856.29 

17651.86 

298.6/028.6 


01-15 


2849.13 

176C7.63 

298.4/028.4 


C 1-30 


2842.02 

17563.70 

298.2/028.2 


01-45 


2834.96 

17520.05 

298.0/028.0 


02-00 


2827.94 

17476.70 

297.8/027.8 


02-15 


2820.98 

17433.64 

297.6/027.6 


02-30 


2814.06 

17390.87 

297.4/027.4 


0 2-45 


2807.19 

17348.41 

297.2/027.2 



21-25N 150-09W 

LONGITUDE WEST 



114-15 

2872.30 

17750.84 

299.5/029.5 

2864.94 
17705.34 

299.3/029.3 

2857.62 
17660. 11 
299.1/029. 1 

2850.35 

17615.17 

298.9/028.9 

2843.13 
17570. 52 
298.7/028.7 

2835.95 
17526.14 

298.5/028.5 

2828.81 

17482.06 
298.3/028. 3 

2821.73 

17438.27 

298.1/028.1 

2814.69 
17394.77 

297.9/027.9 

2807.70 

17351.57 

297.7/027.7 

2800.75 

17308.67 

297.5/027.5 

2793.86 

17266.06 
297.3/027.3 



114-30 

2859.26 

17670.20 

299.6/029.6 

2851.87 

17624.54 

299.4/029.4 

2844.52 

17579.16 

299.2/029.2 

2837.23 

17534.07 

299.0/029.0 

2829.98 

17489.26 

298.8/028.8 

2822.77 

17444.74 

298.6/028.6 

2815.62 

17400.50 

298.4/028.4 

2808.51 

17356.57 
298.2/028.2 

2801.44 

17312.92 

298.0/028.0 

2794.43 

17269.58 
297.8/027. 8 

2787.46 

17226.53 

297.6/027.6 

2780.55 

17183.78 

297.4/027.4 



114-45 

2846.22 

17589.65 
299. 7/029.7 

• 2838.81 
17543.83 
299.5/029.5 

2831.44 

17498.30 

299.3/029.3 

2824. 12 
17453.04 
299. 1/029. 1 

2816.84 

17408.08 

298.9/028.9 

2809.61 

17363.41 

298.7/028.7 

2802.43 

17319.03 

298.5/028.5 

2795.30 

17274.94 

298.3/028.3 

2788.21 

17231.15 

298.1/028.1 

2781.17 

17187.66 
297.9/027.9 

2774.19 

17144.47 

297.7/027.7 

2767.25 

17101.58 

297.5/027.5 
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C THIS PROGRAM COMPUTES TABLES OF POSITIONS 

DIMENSION CPM(36,8)» OBM(36,8) 

READ (5,3) ALATS, ALNGS 
3 FORMAT (2F10.5) 

READ (5,3) ALAT, ALNG 
AMALT = 40. 

ALNG1 = ALNG 
IAB = 34 

7 IF (ALAT. GE. AMALT) GO TO 650 
II = 1 

DO 500 1=1, 12 

DO 400 J = 1, 8 
RALAT = ALAT/57. 29578 
RALATS = ALATS/57. 29578 
IF ( ALNGS. GE. 0.0) GO TO 10 
IF (ALNG. LT. 0.0) GO TO 11 
DLO = ALNGS - ALNG 
ACLO = ABS ( DLO ) 

IF ( ADL0.GE.180. ) DLO = 360. + DLO 
GO TO 13 

10 IF (ALNG. GE. 0.0) GO TO 11 
DLO = ALNGS - ALNG 

IF (DLO. GE. 180. ) DLO = -360. + DLO 
GO TO 13 

11 DLO = ALNGS - ALNG 
GO TO 13 

13 ADLO = ABS (DLO ) 

CONTINUE 

HADLO = C.53ADL0/57. 29578 

TFE = SIN(HADL0)-*SIN(H4DLG)*C0S(RALAT)&C0S(RALATS) 
DL = ALATS - ALAT 
ADL = ABS(DL) 

HADL = 0 . 5’’ : ADL / 57. 29578 

HD = S I N ( H A D L ) ^ S I N ( H A D L ) + THE 

CD = 1. - 2.#HD 

IF (CD.GT.O.O) GO TO 15 

B = -CD 

ARGS = B/SQRT ( 1. - B*-*2) 

RD = 1.5708 + ATAN(ARGS) 

GO TO 16 

15 ARGC = CD/SQRT ( 1 . - CD**2) 

RD = 1.5708 - AT AN ( ARGC ) 

16 CONTINUE 

D = RD * 57.29578 
DNM = 60. # D 
DMS = 6.18 * DNM 
CCL1 = 90. - ABS(ALAT) 

AD = D - CCL1 
DCL1 = ABS ( AD) 

IF ( ALAT. GE. 0.0) GO TO 25 
IF ( ALATS. GE. 0.0) GO TO 28 
GO TO 26 

25 IF ( ALATS. GE.O .0) GO TO 26 
28 CCL2 = 90. + A BS ( A L ATS ) 

GO TO 27 

26 C0L2 = 90. - ABS(ALATS) 

27 CONTINUE 

HC0L2 = C.5$C0L2/57. 29578 
HCCL1 = C . 5 : DCL1/57. 29578 

CA = SIN(HC0L2)*SIN(HC0L2) - S I N ( HDCL 1 ) *S I N ( HDCL1 ) 
CB = 1 . /COS ( RALAT ) 

CC = l./SIN(RD) 

HCCC = CA CB $ CC 
CCCC = 1. - 2. * HCCC 
ARGCC = CCCC/SQRT ( 1 . - CCCC*#2) 

RCCC = 1.5708 - AT AN ( ARGCC ) 

CCC = RCCC * 57.29578 
IF ( ALAT. GE.O. 0) GO TO 30 
IF (DLO.GE.C.O) GO TO 35 
GO TO 36 

30 IF (DLO. GE. 0.0) GO TO 40 
GC TO 37 
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35 CCC = CCC - 180 
GO TO 40 

36 CCC = CCC + 180. 

GO TO 40 

37 CCC = 360. - CCC 
40 CONTINUE 

BCC = CCC + 90. 

IF (BCC. GT. 360 .) BCC = BCC - 360. 

IJ = II + 1 
IK = II + 2 
OPM ( II , J ) = DNM 
OPM ( I J , J ) = CMS 
OPM ( IK, J) = CCC 
OBM ( IK, J ) = BCC 
ALNG = ALNG - 0.25 
4 CO CONTINUE 

ALNG = ALNG1 
II = IK 4 1 
ALAT = ALAT + 0.25 
500 CONTINUE 

WRITE (6,100) 

100 FORMAT ( 1H1 , "NWC",7X, "22.3 KHZ", 2X," 21-49. OS" ,4X, 

1 " 1 14-09 . 8E" ) 

WRITE (6,101) 

101 FORMAT ( "LATITUDE", 47X, "LONGITUDE WEST",47X, 

1 "LAT I TUDE" > 

WRITE (6,102) 

102 FORMAT ( 2 X , "NO RTH" , 5X , " 1 24-0 0" , 7 X , " 1 24- 1 5 " , 7 X , 

1 " 124-30", 7 X, " 1 24-45", 7X," 125-00", 7X , " 125- 1 5" , 7X, 

1 « 125-30", 7X,"1 25-45 ", 9X, "NORTH") 

L = 1 
I CD = 00 

DO 600 K = 1, 36 
GO TO (601,602, 603) L 

6C1 WRITE ( 6,605) I A B , I CD , OP M ( K , 1 ) , 0 P M ( K , 2 ) , OPM ( K , 3 ) , 

10PM(K,4) ,OPM(K ,5) , OPM ( K , 6 ) , OPM ( K , 7 ) , OPM ( K , 8 ) , IAB, ICD 
6C5 FORMAT ( 1H0 , IX , I 2 , "-" , 1 2 , 4X , F8 . 2 , 5X , F 8. 2 , 5X , F 8. 2 , 5X , 
1F8.2,5X,F8.2,5X,F8.2,5X,F8.2,5X,F8.2,8X,I2,"-",I2) 

ICD = ICD + 15 
IF (ICD.EQ.60) GO TO 615 
GO TO 610 
615 ICD = 00 

IAB = IAB + 1 
GO TO 610 

6C2 WRITE ( 6,606) OPM ( K , 1 ) , OPM ( K , 2 ) , OP M ( K , 3 ) , OPM ( K,4 ) , 

1 OPM ( K , 5 ) , OPM ( K , 6 ) , 0 PM ( K , 7 ) , OPM ( K , 8 ) 

606 FORMAT ( 1 IX , 8( F 8 . 2 , 5X ) ) 

GG TO 610 

603 WRITE (6,607) OP M ( K , 1 ) , OBM ( K , 1 ) , OP M ( K , 2 ) , OBM ( K , 2 ) , 

10PM(K, 3) ,OBM(K ,3) , 0PM(K,4) ,0BM(K,4) ,OPM(K ,5) ,0BM(K,5) , 
1 OPM ( K , 6 ) , 03M(K ,6) , OPM ( K , 7 ) , OBM ( K , 7 ) ,GPM< K ,8) ,0BM( K,8) 
6C7 FORMAT ( 9X , 8 (F 5. 1 , "/" , F5 . 1 , 2X ) ) 

GC TO 610 
610 L = L + 1 

IF (L.GT.3) L = L/3 
600 CONTINUE 
GO TO 7 

650 WRITE (6,651) 

651 FORMAT (1H1) 

END 
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